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Q^ . ABSTRACT 

0^ \ A theory is developed for the dynamics of eccentric perturbations [oc exp(±i0)] of a disk galaxy residing 

in a spherical dark matter halo and including a spherical bulge component. The disk is represented as a 

P5 ' large number N of rings with shifted centers anrf with perturbed azimuthal matter distributions. Account 

^ \ is taken of the dynamics of the shift of the matter at the galaxy's center which may include a massive 

black hole. The gravitational interactions between the rings and between the rings and the center is 

Q^ ' fully accounted for, but the halo and bulge components are treated as passive gravitational field sources. 

Equations of motion and a Lagrangian are derived for the ring+center system, and these lead to total 

energy and total angular momentum constants of the motion. 

^ \ We first study the eccentric motion of a disk consisting of two rings of different radii but equal mass, 

J^ . Md/2. For small Md the two rings are stable, but for Md larger than a threshold value the rings are 

unstable with a dynamical time-scale growth. For Md sufficiently above this threshold, the instability 

acts to decrease the angular momentum of the inner ring, while increasing that of the outer ring. The 

pi. • instability results from the merging positive and negative energy modes with increasing Md- Secondly, 

Qs I we analyze the eccentric motion of one ring interacting with a radially shifted central mass. In this case 

O^ . instability sets in above a threshold value of the central mass (for a fixed ring mass), and it acts to 

increase the angular momentum of the central mass (which therefore rotates in the direction of the disk 

matter) , while decreasing the angular momentum of the ring. 

Thirdly, we study the eccentric dynamics of a disk with an exponential surface density distribution 

O ' represented by a large number of rings. The inner part of the disk is found to be strongly unstable. 

Angular momentum of the rings is transferred outward and to the central mass if present, and a trailing 






Oh' 



U 



C/3 ■ one-armed spiral wave is formed in the disk. Fourthly, we analyze a disk with a modified exponential 



a 



X 



density distribution where the density of the inner part of the disk is reduced. In this case we find much 
slower, linear growth of the eccentric motion. A trailing one-armed spiral wave forms in the disk and 
becomes more tightly wrapped as time increases. The motion of the central mass if present is small 
compared with that of the disk. 

L INTRODUCTION Peak 0.9 meter observations of NGC 1073 (Kornreich et 

-V ,,, , , J. r ■ 1 1 • 1 al. 1998) indicate that even the optical center of light 

Although studies of spiral galaxies commonly assume , , • , i r , , , r i • i • , , 

, . .,.1 • ? . -11 i u J u may be displaced from the center oi the optical isophotes 

an axisymmetric equilibrium state, possibly perturbed by ^ f. . . r ,■, , • . .„ , „ f. „ 

. , J.U • ■ -J j-u J- 1 ol the maior part of the galaxy m a sigmncant traction OI 

spiral arms, there is growing evidence that many galax- j f t. j t> 

ies lack this symmetry. Based on optical appearance, ap- „. ' . . r;, , . , • ,. , , . 

i 1 onoy r J- 1 1 • 1 -i^-i • -c i- ui i he origin ot the lopsideness m disk galaxies remains un- 

proximately oO/o oi disk galaxies exhibit signmcant lop- , . tV , i • , 1 /1r^o/^^ i ■ i i • 

., 1 )) /o- c ry -,1 innc; t^ • u tr i> Certain. Baldwin ei a/. (1980) proposed a simple kinematic 

sidedness (Rix & Zaritsky 1995; Kornreich, Haynes, & ... ,.,,.„. ) ■ r ,^ ^ , 

T 1 -inno\ ^ ■ I. j- j-i, i c j- -f o u model m which ditterent rmgs ot the galaxy, assumed non- 

Lovelace 1998) which supports the early hndmgs ot Bald- . , ,. • •;• „ , • n , n ;, • . i 

T J A 11 c o • • rinan\ \ rn(» interacting, are initially shitted irom their centered equi- 

wm, Lynden-Bell, & Sancisi (1980). As many as ~ 50% ,., . °'.,. mi i -n i • • .i n 

r • 1 1 • 1 1 , r .u i J librium positions. Ihe shitted rmgs precess m the overall 

ol spiral galaxies show departures trom the expected sym- •, ,r , , ;• , . ,. ,. .; ; ,, , r 

, • , , J 1 u 1 TJT T ci /ID- u-t PC gravitational potential m a direction opposite to that ot 

metric two-horned global HI line pronle (Richter & San- ° , • t-. ,i ■ , i 

• • inn/f tr 7 i 1nno^ t7 iu ttt c the mass motion. Because the precession rate decreases m 

cisi 1994; Haynes et al. 1998). l<urthermore, HI maps of , .,, ,•,,•; ^ ..,.,,., , , , 

1 , • r 1 ATnri oeoi n^ 1nn'7^ general with radial distance, an initial disturbance tends 

several galaxies, tor example, NCL 3631 (Knapen 1997), f « ■ i „ ■ . , i- • , 

NGC 5474 (Rownd, Dickey, & Helou 1994), and NGC 7217 *° J'f ^P ^^*° ^J'^urf- '^''t^,"^ '"".^ 1 "iT.t 

(Buta et al. 1995), have revealed offsets between the op- '^^^ ^^'' *^^", *^^ ^^^^^^ *™^- ^'^^^' ^f^ .^™*^ [^^^^ 

tical centers of lie:ht and their kinematic centers Two ^^^^> ^^^^ ™^*^'^ extensive computer simulation studies of 

, r, . ,• 11 1 -1 J 1-1. i.1 the unstable eccentric motion ot matter m the nuclei ot 

examples ot kmematically lopsided galaxies have recently , . , . , , , , . , . , , , n- 

1 ,. J u c J- \ 1 nc\no\ n.i.1, 4. u galaxies which they suggest is pertinent to the on-center 

be discussed by Swaters ei a/. (1998). Other recent obser- , • , ,. i ^ ^ ■ r i uto^ 

,. 1 ii rr tLj r> rri, i ,• iiuclci obscrvcd lu a numbcr ot galaxics, tor cxample, JV131, 

vations such as the Bubble Space lelescope observations _ _„„ i ^1^1 jr- i j 

M33, and MlOl. 



of the nucleus of M31 (Lauer et al. 1993) and the Kitt 



Understanding of the origin of the observed disk asym- 



metries is important because it provides clues to ongoing 
accretion of gas and the distribution of dark, unseen mass 
in the halo. Schocnmakers, Franx, and de Zeeuw (1997), 
interpret optical asymmetry as an indicator of asymme- 
try in the overall galactic potential, and therefore an in- 
dicator of the spatial distribution of the dark matter in a 
galaxy which may have a triaxial distribution. By analyz- 
ing the spiral components present in the surface brightness 
or H I distribution, the velocity gradient, and therefore, 
the shape of the gravitational potential, may be uncov- 
ered. Jog (1997) has studied of the orbits of stars and gas 
in a lopsided potential, and shows that lopsided potentials 
arising from disks alone are not self-consistent; rather, a 
stationary lopsided disk may be responding to asymme- 
tries in the halo. 

Zaritsky and Rix (1997) proposed that optical lopsided- 
ness arises from tidal interactions and/or minor mergers. 
Such mergers arc often suggested as the most likely con- 
tributors to galaxy asymmetry, even when no interacting 
companions are evident. While galaxies such as NGC 5474 
are well-known to be under the tidal influence of neighbors, 
the apparently long-lived kinematic offsets of other rela- 
tively isolated objects, and the common asymmetries in 
flocculent (as opposed to tidally induced "grand design") 
spiral galaxies, are not explained by simple tidal inter- 
action models, which produce only transient asymmetric 
features. 

A further possibility is that an optical disk may be in a 
quasi-stationary lopsided state in a symmetric potential, 
as discussed by Syer and Tremaine (1996). In this model, 
gaseous and stellar matter swirl about the minimum of 
the halo potential in a state not fully relaxed. The result 
is a lopsided flow within a symmetric mass distribution. 
Numerical simulations of this situation have been done by 
Levine and Sparke (1998) using a gravitational A^— body 
tree-code method (see Barnes & Hut 1989) for disk galax- 
ies shifted from the center of the main halo potential. The 
results are suggestive of lopsidedness with large lifetimes. 
An A^ — body simulation study of a rotating spheroidal stel- 
lar system including the dynamics of a massive central ob- 
ject by Taga and lye (1998a) indicates that the central 
object goes into a long lasting oscillation similar to those 
found earlier by Miller and Smith (1988, 1992) and which 
may explain asymmetric structures observed in M31 and 
NGC 4486B. A linear stability analysis of a self-gravitating 
fluid disk including a massive central object also by Taga 
and lye (1998b) indicates a linear instability (Taga & lye 
1998b). We comment on the relation of this work to the 
present study in the conclusions section of this work. 

Here, we develop a theory of the dynamics of eccentric 
perturbations of a disk galaxy residing in a spherical dark 
matter halo. We represent the disk as a large number A^ of 
rings as suggested by Baldwin et al. (1980) (and Lovelace 
(1998) for the treatment of disk warping). In contrast with 
Baldwin et al., the gravitational interactions between the 
rings is fully accounted for. We show that for general ec- 
centric perturbations, the centers of the rings are shifted 
and the azimuthal distribution of matter in the rings is 
perturbed. The ring representation is analogous to the 
approach of Contopoulos and Gr0sbal (1986, 1988), where 
self-consistent galaxy models are constructed from a finite 
set of stellar orbits. 



Section 2 develops a theory for treating eccentric per- 
turbations of a disk galaxy. The assumed equilibrium is 
first discussed (§2.1), and a description of the disk pertur- 
bations is developed (§2.2). The representation of the disk 
in terms of a finite number A^ of rings is presented (§2.3), 
and the ring equations of motion are derived (§2.4). We 
renormalize the ring equations so as to reduce the nearest 
ring interactions (§2.5). The dynamics and influence of 
the displacement of the center of the galaxy, which may 
include a massive black hole, is discussed separately (§2.6). 
We obtain an energy constant of the motion for the dynam- 
ical equations (§2.7), the Lagrangian, and the conserved 
total canonical angular momentum (§2.8). 

We discuss the nature of the precession of a single ring 
in §3. In §4 we study the eccentric motion of a disk con- 
sisting of two rings and show that this motion is unstable 
for sufficiently large ring masses. In §5 we study the ec- 
centric motion of one ring including the radial shift of the 
central mass and show that this situation is unstable for 
sufficiently large mass of the center and/or of the ring. 
Section 6 presents numerical results for the eccentric dy- 
namics of disk of many rings including the radial shift of 
the central mass. Section 7 summarizes the conclusions of 
this work. 

2. THEORY 

2.1. Equilibrium 

The equilibrium galaxy is assumed to be axisymmet- 
ric and to consist of a thin disk of stars and gas and a 
spheroidal distributions consisting of a bulge component 
and a halo of dark matter. We use an inertial cylindrical 
(r, (/), z) and Cartesian (x, y, z) coordinate systems with the 
disk and halo equatorial planes in the z — plane. The 
total gravitational potential is written as 



$(r, z) = $d + $6 



$/, 



(1) 



where $d is the potential due to the disk, $;, is due to the 
bulge, and ^h is that for the halo. The galaxy may have 
a central massive blac k hole of mass M;,^ in which case a 
term $(,/i = —GMiih/^r"^ + z^ is added to the right-hand 
side of (1). The particle orbits in the equilibrium disk are 
approximately circular with angular rotation rate ri(r), 
where 

^^* =nl + nl + nl. (2) 

The equilibrium disk velocity is v = rQ(r)(p ■ A central 
black hole is accounted for by adding the term 17^^ = 
GMbh/r^ to the right-hand side of (2). 

The surface mass density of the (optical) disk is taken 
to be Sfi = E(joexp(— r/rc() with S^o and r^ constants and 
Md — 2TTr'^T^do the total disk mass. The potential due to 
this disk matter is 



r or 



$d(r,0) = - 



GMd 

Td 



R[Io{R)Ki{R) - h(R)Ko(R)] 



and the corresponding angular velocity is 



^2 



l^[I,{R)Ko{R}-h{R)K,{R)] , (3) 



where R = r/{2rd) and the I's and K's are the usual mod- 
ified Bessel functions (Freeman 1970; Binney & Tremaine 



Fig. 1. — Sample disk rotation curve v^{r) for the values M^ = 6 x IO'^'^Mq and r^ — 4 kpc for the disk, M^ = 5 x IO^Mq 
and ri, — I kpc for the bulge, and Vh = 250 kni/s and r/i = 5 kpc for the halo, using expressions given in §2.1. 



1987, p.77). Typical values are M^ = 6 x IO^OMq and 
rd — 4: kpc. For these values, Vd ^ ^jGMd/rd ~ 255 
km/s. 

The potential due to the bulge component is taken as a 
Plummer model 



$h = - 



GMh 



{rl+r^ + z-^y/-^ ' 



where Mf, is the mass of the bulge and r;, is its character- 
istic radius (Binney & Tremaine 1987, p. 42). We have 



o2 ^ GM, 

' (r2+r2)3/2 



(4) 



Typical values are M^ — 10 M© and r^ — 1 kpc, and for 
these values v\j = ^jGM\Jr\, ~ 208 km/s. 
The potential of the halo is taken to be 



$/, 



A ln(r^ + r^ + z") , 



where v^ = const is the circular velocity at large distances 
and Th = const is the core radius of the halo. We have 



ni 



r? + r^ 



(5) 



Typical values are v^ ^ 200 — 300 km/s and r}^ ^ 2 — 20 
kpc. Figure 1 shows an illustrative rotation curve. 

2.2. Perturbations 

We treat the disk as fluid and use a Lagrangian rep- 
resentation for the perturbation as developed by Frieman 
and Rotenberg (1960). The position vector r of a fluid 
element which at t = was at ro is given by 



r = ro +^(ro,i) 



(6) 



That is, ro(i) is the unperturbed and r(t) the perturbed 
orbit of a fluid element. This description is applicable 



to both the disk gas and the disk stars which are in ap- 
proximately laminar motion with circular orbits. The per- 
turbations of the halo and bulge are assumed negligible 
compared with that of the disk. For these approximately 
spheroidal components, the particle motion is highly non- 
laminar with criss-crossing orbits with the result that their 
response "averages out" the disk perturbation. 

Further, the perturbations are assumed to consist of 
small in-plane displacements or shifts of the disk matter. 



i^irY + ^4,<^ 



(7) 



with azimuthal mode number m — 1. That is, S,r and ^^ 
have (?!)-dependences proportional to exp(«0). 

From equation (6), we have v(r, t) = vo(ro, t) + d^/dt + 
(v • V)$. The Eulerian velocity perturbation is 5v(r,t) = 
v(r,t) — vo(r,t). Therefore, 



(5v(r, t) 



dt 



+ (vV)^-(^-V)v 



The components of this equation are 



(8) 



(9) 



where 



-4-"<4 



and n' = dn/dr . 

The main equation of motion is 



dSv 
dt 



<5F 



-V(5$ 



(10) 



where SF is the perturbation in the gravitational force (per 
unit mass), and 6^ is the perturbation of the gravitational 
potential. The pressure force contribution is small com- 
pared to SF by a factor {vth/v^,)"^ ^ 1 and is neglected, 
where Vth is the 'thermal' spread of the velocities of the 
disk matter. Also, note that 



d 5-v d Sv 



dt 



dt 



+ (v V) (5v + ((5v • V) V 



(11) 



The components of (11) give 



dt 

dSv 
dt 



= VSv^ + {Kl/2n)SVr , 



(12) 



where k^ = (l/r^)d{r'^il'^)/dr is the radial epicyclic fre- 
quency (squared). 

The perturbation of the surface mass density of the disk 
obeys 






-V • (S (5v + SY. v) 



where S(r) is the surface density of the equiUbrium disk. 
Because V • (S vq) = 0, this equation impUes 

(5S = -V • (S $,) . (13) 

The perturbation of the gravitational potential is given by 

, (5E(r',i) 



(5$(r,t) = -G f d^r' 



(14) 



where the integration is over the surface area of the disk. 

2.3. Ring Representation 

We represent the disk by a finite number TV of radi- 
ally shifted plane circular rings. The matter distribution 
around each ring is also perturbed. An elliptical distor- 
tion of a ring corresponds to an to = ±2 which is not 
considered here. This description is general for small shifts 
(dr/r)'^ ^ 1, where the linearized equations are applicable 
(Lovelace 1998). Of course, the orbit of a single perturbed 
particle is not in general closed in the inertial frame used. 
However, the orbit is closed in an appropriately rotating 
frame, and this rotation rate is simply the angular preces- 
sion frequency of the ring lo discussed below (Baldwin et 
al. 1980). The disk is assumed geometrically thin. 

For the equilibrium disk we take 



S(r) 



M, 



N 

E 



. , 27rr\/ 27rAri 



exp 



(r- 



(15) 



where Mj is the mass of the j ring, rj is its radius with 
< ri < r2 ... rjv, and Ar^ ^ rj is its width. The motion 
of the central part of the disk (r < ri ) is treated separately 
in §2.6 

A physical choice for the rings distribution will have the 
ring spacing of the order of the disk thickness, rj+i — rj — 



0{lS.z). Further, we assume [vj+i 



<C Tjj^irj m or- 



der to simplify the calculation of the ring interaction as 
discussed in Appendix A. For example, a possible choice 
is Tj = ri + {j — l)Sr with ri = 1 kpc, Sr = 0.5 kpc, and 
Mj = 2TTrj5rT,dirj). For Avj = Sr /V2h[2 « Sr/ 1.177, 
the profile of a ring falls to half its maximum value at Sr 
so that equation (15) gives a fairly smooth representation 
of Mr)- 



The perturbation in the disk's surface density is 

5S(r, 0, t) - ^ <5S, , (5S, = -V • (S,C) , (16) 



where Ej = (Mj/27rr\/27rArj) exp[-(r - rj)'^/{2Ar^)]. 

We can express different moments of the perturbed disk 
in terms of the rings. For example, the center of mass of 
the disk is 

<r>=^^</^>, (17) 



EM, 



where 



Mj < r, > 



d'^x r (5Eo 



- ^''JT.^> 



(18) 



where an integration by parts has been made. 
We can write in general 



O0 



^OV ' 



Sjx sin (h + Sjy cos 



ov ' 



(19) 



Here, ejx,y and Sjx,y are the ring displacement amplitudes: 
ejx.y represents the shift of the ring's center, and Sjx,y rep- 
resents in general both the shift of the ring's center and 
the azimuthal displacement of the ring matter. Firstly, 
notice that for ejx,y = and Sjx,y ^ 0, there is no shift of 
the ring's center but rather an azimuthal displacement of 
the ring matter. In this case (5Sj = —{'Sj/rj){d£,j^/d(j)). 
Secondly, notice that Sjx,y = ^jx,y corresponds to a rigid 
shift of the ring without azimuthal displacement of the 
ring matter. For example, a rigid shift in the x— direction 



nas ejx — ^jx anci ejfj 



= Sjy SO that S^j 



and S^j 



-ejxsincj). In this case, V-^„ = so that 



(5Sj = -V-(Sj^j) = -Crj{dT,j/dr). Figure 2 shows the 
nature of perturbations with a shift of the ring's center 
and with an azimuthal displacement of the ring matter. 
Equation (18) now gives 






(20) 



For the case of a rigid shift of a ring, ejx^y = Sjx^y, the 
center of mass position is simply < rjx.y >~ £jx,y, as ex- 
pected. 

Similarly, the velocity perturbation of the disk can be 
written as 

EM,<^v,>^ (21) 



where 



< Sv >-- 



1 



EM, 



"^ '^^^ ^" M ' ^'^^■'^ 



Sj(5vj) 



(22) 



Evaluation of (22) gives 

< Svjx > 



< 5Vjy > 



' Ojx 



'jy + "jy 



(23) 



The influence of the eccentric motion on the rotation 
curves is discussed at the end of §6.3. 



Fig. 2. — Drawing of two perturbed rings with equilibrium radii ri — 1 and r2 = 2 in arbitrary units. The center of 
ring 1 is rigidly shifted in the x— direction from the origin by eix — 0.4 and eiy = 0. The center of ring 2 is at the origin 
while the distribution of matter around the ring is perturbed as indicated by the small circles. For this ring, S2x = 0.4 
and (52 J, = 0. 



2.4. Ring Equations of Motion 

Equations (12) and (19) give the ring equations of mo- 
tion, 

'ix + 2ney ~ n^Ex - 2n{5y - nSx) = < sf^ > , 

■iy - 2^6,; - n^ey + 2fi(4 + ilSy) = < S F^^ > , 

4 + 2^4 - n'^Sx - 2n{ey - flex) = - < SF^ > , 



Sy - 2nSx - n^6y + 2n{ex + ney) = < &f9 > 



(24) 



where the j subscripts are implicit, where the angular 
brackets indicate the average over the ring < (..) >= 
27r/ rdr(..)Sj(r)/Mj, where Cl^ = fl^ - 2nrn' , and 



SF, 



c.s 



[cos (j), sin (p] 5 Fa , 



with a = r, (j). 

We now evaluate the different force terms on the right- 
hand side of (24). For this it is useful to write (5S = 
SY^a + <5Sfc, where 

oi^a = ^ , oSfc == HI—' (25) 

r or r 0(p 

from equation (13). The corresponding contributions to 
the potential, 5^ = S^a + S^b evaluated at (r, (p) are from 
(13), 



(5$Q (r, (/)) = G ^ Mfc {ekx cos (j) + euy sin < 



r' dr' 



r'S{r'\rk) 



d^' cos * 



27r R{r,r') 



(26) 



and 



(5$b(r, 0) = -G^MkiSkxCos<l) + Sky sin (j)) x 



' , , , S{r'\rk) rd^ cos* 
r dr '*• 



10 r' J 27r R{r,r') 

where R^{r,r') = r^ + (r')^ - 2rr'cosvI', 



(27) 



S{r\rk) = 



_ 27rSfe(r) _ cxp[-(r-rfe)V2Ar 



Mfc 



r\/27rArfc 



and * = 0' - (/). 

Evaluation of the force on the j*"^ ring due to the other 
rings gives 

Mj < 5F^j > = XI ^'^ok'^kx + DjkSkx) , 

k 

Mj < 6F^j > ^ Yl (^3k£ky + DjkSky) , 
k 

-Mj < SF^^ > ^ Y. i^ikSkx + D'^^ekx) , 
fc 

M,<5Fg> = X(ii;,fc4,+i?;fcefe,), (28) 

k 

where the 'tidal coefRcients' are 

C.k = -GMjMk [ [ rdr r'dr' x 



S{r\r,) 



d[r'Sir'\rk)] a/C(r,r') 



dr' 



(29) 



D,k = GM.Mk / / rdr r'dr' x 



IjlVlk 



5(r|r,)5(r |rfc)— j^— 



(30) 



D'j^ = ^GMjMk / / rdr r'dr' x 
5(r|rj)5[r'^(r'|rfe)] 



'dr' 



Kfe = GMjMk / / rdr r'dr' x 



/C(r,r'), 



(31) 



Sir\r,)Sir'\n) ^^ 



where 



/C(r,r') = 



d* cos* 

2^i?(r,r') ' 



(32) 
(33) 



and where the r, r' integrals are all from to oo. 
Formal integration by parts of (29) gives 

C,k = GM.Mk / / rdr r'dr' x 



-■]k 



■IjIVlk 



5(r|r,)^(r'|rfe) 



d^lC{r,r') 
dr' dr 



(34) 



which shows that Cjk ~ Gkj- Also, integration by parts of 
(31) gives 



D^fe = GMjMk I I rdr r'dr' x 



5(r|r,)^(r'|r,) 



dlC{r,r') 
r dr' 



(35) 



and this shows that Djk = -DL- Expressions for the tidal 
coefficients in terms of elliptic integrals are given in the 
Appendix. 

Following the approach of Lovelace (1998), we introduce 
the complex displacement amplitudes 

£j = Ej^ - iejy ^ ej{t) cxp[-itpj{t)] , (36) 

Aj = 6jx - i6yy = 6^{t) exp[-z'(/'j(t)] . (37) 

Here, e^ > is the amplitude of the shift of the ring's cen- 
ter, and ipj is angle of the shift with respect to the x— axis; 
Sj > is the amplitude of the azimuthal displacement of 
the ring matter, and ipj is the angle of the maximum of 
the ring density also with respect to the a;— axis. If ip{t) 
and tpit) increase with time, the ring precesses in the same 
sense as the particle motion and we refer to this as forward 
precession. The opposite case, with ip{t) and ip(t) decreas- 
ing with time, is termed backward precession. 



We now combine equations (24) and (28) to obtain the 
ring equations of motion, 

Mj ( £j + 2inj£j - n'^jSj - 2injAj + 2ri|Aj') = 

J2iCjkSk + DjkAk) , (38) 



Mj (Aj + 2injAj - n'^jAj - 2inj£j + 2n'^j£j ] = 



7 ^ [EjkAk + D'jf.£k 



(39) 



where j — \..N , Vtj = ri(rj), and Vtj = ^(r^). 

2.5. Renormalization of Ring Equations 

Here, we redo the ring equations of motion so as to di- 
minish the strong tidal interactions of nearest neighbor 
rings due to the terms ex Gj,j+k- First, we rewrite the 
right-hand side of (38) as 

Y,[C,k{£k - £,) + D,k{Ak- £,)] + 

k 

£iY.^C,k+D,u), (40) 

k 

Similarly, we rewrite the right-hand side of (39) as 
^ [E,k (Afe - A,) + D'^^ {£k - A,)] + 



^E(^: 



jk 



D'.: 



Jk) 



We define 



02 



ni 



jf X! (^j*= + -^-j'*^) 



(41) 



(42) 



Using the relations fi^ = fi^ _ r{dQF-/dr), fi^ = fi^ + fig 
fi^, and equation (A14) of the Appendix gives 

fi^^. = fi^fr- 

(If there is a central black hole Mfy^ then the right- 
hand side of (43) also has the term ~[r{d^1f^/dr)\j — 
SGMbh/r^.) Similarly, we define 



Q-r +jr2^(^Jk + E,k). (43) 



n 



Aj 



ni 



^X^D'^k+E,,). 



(44) 



The ring equations of motion now become 



(j ( £j + 2iVLj£j - 


- Vl%£, - 2zr!, A, + 2Vl]A, 


y, [Gjk {£k - 

k 


-£^)^D,k{A^-£,)\ , 


'j (Aj + 2injAj - 


- fiij Aj - 2ifij4- + 2f^jfj 



^[ii;,,(Afc-A,) + i?;,(ffe-A,)] 



(45) 



(46) 



In the large N limit, the sums in equations (45) and (46) 
go over to bounded integrals. The diagonal elements, Cjj 
(and -Ejj), are absent from (45) and (46). Note also that 
the self-interaction of a ring vanishes as it should for the 
case of a rigid shift where £j — Aj . 



2.6. Dynamics and Influence of "Center" 

The dynamical equations (45) and (46) do not account 
for the part of the disk inside the innermost ring ri. Also, 
there may be a massive black hole M^h near the galaxy 
center. We treat this central region separately as a point 
mass Mq, 



Mo = Mbh + 27r / rdr Sd(r) , (47) 

Jo 

where r'l = ri — Sri/2. 

The horizontal displacement of the "center" is 

eo(i) = eoxX + eoj/Y • (48) 

The equation of motion for cq taking into account the ec- 
centric displacemnts of the rings is 



Mo 



where 



d^£. 



N 



1 



—^ + ^IqEq j — - 2^ Fok ( £k - xAfe 



£o = eoa; - ieoy = eo cxp{-iipo) 



(49) 



(50) 



is the complex displacement amplitude of the "center" ; 
Fok = GMoMk/rl, k = l..iV, are the tidal coefficients 
between the "center" and the rings; and 



ni^[Y-l\ ^('J^'^^ + '^^) 







r 



dr 



2M, 



h:T.Po, (51) 



is the angular oscillation frequency of a particle at the 
galaxy center. From §2.1, we have [(l/r)(9$b/i9r)]o = 
GMb/rl and [{l/r){d<^h/dr)]o = vl/rl- For all of the 
considered conditions, we find ri§ > 0. 

The influence of the displaced "center" on the rings is 
included by adding the force terms, due to the "center's" 
displacement, to the right-hand sides of equations (45) and 
(46), 



Mi £. 



— — 2i^0j"^0 + 



(52) 



Mj A, 



= + FojSo + ■■ , (53) 

where the ellipses denote terms in equations (45) and (46). 



-Ejk{Sj 



lMo{el + nlel) + J2FoJe„-e,-leo-S,) , (54) 



and where the real vectors e ~ e^^+Syy and d = Sxi+Syy 
are useful here. 



2.8. Lagrangian 

By inspection, we find the Lagrangian for the ring sys- 
tem C{ejx,(^]x,--), 

^ = i E ^^^ Is' + ^h^' + ^' + ^h^' ^^h ■ ^^ 

-2il,[lej - Si) X (ej - S/)] ■ i] + 1a/„(«; - ill4) 



-IT.T. 






Sk? 



Cjki^j — ^k) +'2Djk{ej — Sk) 



Because dC/dt = 0, the Hamiltonian 

dC 



^-E^ 



dijoo 



C 



(55) 



(56) 



is a constant of the motion. It is readily verified that 
H = E. 

We can make a canonical transformation {ejx,ejy) — > 
{ej,ipj), (SjxjSjy) — > {Sj,ipj) to obtain the Lagrangian 
as C — C{ej,Sj,ej,6j,ipj,ipj). Note for example that 
e — > e^ + e'^ifi'j. It is then clear from the azimuthal symme- 
try of the equilibrium that C is invariant under the simul- 
taneous changes cpj -^ (pj + 9, ipj -^ ijjj + for j/ = 1, .., iV, 
where 9 is an arbitrary angle. Thus 



E 



dL dC 



= 



and consequently the total canonical angular momentum 
of the ring system. 



2.7. Energy Conservation 

An energy constant of the motion of the ring system 
can be obtained by multiplying equation (52) by £* and 

(53) by A* (with [...]* denoting the complex conjugate), 
adding the two equations, summing over j, and dividing 
by 2. (This factor of 2 makes the kinetic energy of a rigidly 
shifted ring with Ej = Sj equal to (l/2)Afje|.) Further, 
we multiply equation (49) by £q and add the result to the 
previous sum. In this way we find dE/dt = 0, where 



E- 



jE«* 






^,^ 



m'^e.-S, 






ky^j j 



^jki^j — ^k) + '^Djk{ej — 8k 






(57) 



is another constant of the motion. Evaluating (57) gives 



\E», 



e^Jifi, - rtj) + S^Jijj - flj) 



jWj 



+ 2QjejSj cos{ipj — Tpj) 



Moe„(fo 



(58) 



where the last term represents the angular momentum of 
the galaxy center. The last term within the square brack- 
ets can also be written as 2D,j{ej ■ Sj). 

The constants of the motion Ti and V^ are valuable for 
checking numerical integrations of the equations of motion 
(52) and (53). 



Fig. 3. — The top panel (a) of the figure shows the radial dependence of the four frequencies Wq (a = 1,..,4) of the 
modes of oscillation of an eccentric ring normalized by the disk's angular rotation rate i^(r). The bottom panel (b) shows 
the corresponding radial dependence of the mode amplitude ratio 5/A as discussed in the text. The case shown is for 



the galaxy parameters of Figure 1 with innermost ring of radius ri 
using equations (All), (43), and (44). 



1 kpc. The values ft 



J' 



rjgj, and n^j 



were obtained 



3. ECCENTRIC MOTION OF A SINGLE RING 

Consider the eccentric motion of a particular ring with 
the other rings not excited. This is not a self-consistent 
limit because gravitational interactions will in general ex- 
cite all of the rings. However this limit is informative. 
With £j ex Aj ex exp(— iwi), j = 1, .., A^, where uj the ring 
precession frequency, we get 



[{uj -fl)^ + E]£ + [2n{uj -n)+d]A = , 

[2n{Lu -n) + d]£+[{uj-n)^ + D]A = o , (59) 

where the j subscripts arc implicit, E = ill, — fJ? 



'£] 



Dj,/M„ D = ni^ ~ nj - D,,/M„ and d = D„/M,. 
With w; = cj — ri, we get 



{D + E - AQ' 



)w^ - 4dQw 



ED-d^ = 



(60) 



For the limit of many rings, the terms involving d become 
negligible, and (60) can be readily solved to give four real 
roots iiO < ED < (An^ -D~ EfjA. 

Figure 3a shows the radial dependence of the four ring 
precession frequencies Wq, {a = 1,..,4) corresponding to 
the four modes of oscillation. These modes are the ana- 
logues of the normal modes of vibration of a non-rotating 
system (see Lovelace 1998). 

Three of the modes in Figure 3a have positive frequen- 
cies (wq > for a = 2,3,4) so that they have forward 
precession with 



•fa 



LUat + const , and -tpa = ^at + const' . (61) 



These relations follow from (36) and (37) because |£| = e 
and |A| = (5 are constants. The other mode [a — \) has 
a;i < and therefore backward precession. 



Figure 3b shows that the two modes a = 2, 3 have 8/ A 
small compared with unity near the inner radius of the 
disk, ri. Note that £/A — corresponds to a pure az- 
imuthal shift of the ring matter without a shift of the ring 
center. The mode a — 1 with backward precession has 
£ I A ^ 1. As mentioned, £ / A = 1 corresponds to a rigid 
shift of the ring center. 

Evaluation of the ring energy for the four modes using 
equation (54) shows that the modes a = 1, 2 have negative 
energy whereas the modes a = 3, 4 have positive energy. 
The negative energy modes are unstable in the presence of 
dissipation, for example, the force due to dynamical fric- 
tion (sec for example Lovelace 1998). 

Note that for vanishing ring mass (D ex d ex Mj -^ 0), 
the middle two roots approach w = fl±[DE{'iQ'^ - E)]''-^'^ . 
Thus for Mj -^ 0, there are only three different roots of 
equation (60), uj ^ n and uj ^ n± (4^7^ _ S)i/2_ 

Outside of the central region of a galaxy we have fl ~ 
1/r. Consequently, the radial dependences of the mode 
frequencies Uair) ex ± 1/r will tend to "wrap up" an ini- 
tially coherent asymmetry into a tightly wrapped spiral 
(in the absence of ring interactions). The forward pre- 
cessing modes (a = 2 — 4) will give a trailing spiral wave, 
(fa cx; 1/r, (with respect to the azimuthal motion v^ > 0), 
whereas the backward precessing mode (a = 1) will give a 
leading spiral wave, (pi ex —1/r. The case of mode a = 1 
for non-interacting rings was discussed earlier by Baldwin 
et al. (1980). 



4. ECCENTRIC MOTION OF A "DISK" OF TWO RINGS 

Here, we consider the eccentric motion of a "disk" con- 
sisting of two interacting rings, one of mass Mi and radius 
ri, and the other of mass M2, radius r2. The values of 



^2 



ri\- {j — 1,2) are given by equations (All), 



(43), and (44) with the bulge and halo potentials as given 
in §2.1. Thus the present treatment is self-consistent (in 
contrast with the previous subsection). With £j and Aj 
proportional to exp(— zwi), equations (45) and (46) give 

[ - (t^ - nif + nl- f7|i + Dn]£i 

+ [2ni{ni -uj)-Dn]Ai 
= [C'i2{£2 -£i)+ Du{A2 - fi)]/Mi , (62) 

[-(oj- ni)^ + nl- nil + ^n] Ai 

+ [2ni{ni - Lu) - D[i]£i 

= [Ei2{A2 - Ai) + D[^i£2 - Ai)]/Mi , (63) 



[-{uj-n 



'^2 ^^£2 



D22]£2 

+ [2r!2(r!2-tj)-D22]A2 
= [Ci2(£:i -£2) + ^'12 (Ai - A2)]/M2 , 

[-{u-n2f + nl-n\^ + D'^^]A2 

+ [2^2(^2 - uj)£2 - D'^^] 

^ [^12 (Ai - A2)+Di2{£l - A2)]/M2 . 



(64) 



(65) 

For a non-zero solution, the determinant of the 4x4 ma- 
trix multiplying (£1 , Ai, £2, A2) must be zero. This leads 
to an eighth order polynomial in uj which can be readily 
solved (with Maple R. 5) for the frequencies of the 8 modes 
{a^ 1..8). 

Figure 4 shows the behavior, including instability, of 
a system of two rings of equal mass Mi = M2 so that 
the "disk" mass is Md = 2Mi. For Md ^ 0, the modes 
a = 3,4 approach ^,2 and the modes a = 5,6 approach 
r^i which agrees with the behavior found in §3. Note that 
for small Md modes a = 1 — 4 are associated with ring 2 
while modes a = 5 — 8 are associated with ring 1. In the 
absence of interactions, the rings are stable independent of 
their mass. As the disk mass increases, the modes a = 4, 5 
approach each other and merge at Md ~ 0.276 x lO^^M© 
and give instability with Xm{uj) = Wi > as shown in Fig- 
ure 4b. We refer to this as the "first instability." Notice 
that the onset of instability corresponds to a merging of the 
positive energy mode a = 4 of ring 2 with the negative en- 
ergy mode a = 5 of ring 1 (see §3). The interaction of pos- 
itive and negative energy modes is a well-known instabil- 
ity mechanism (see for example Lovelace, Jore, & Haynes 
1997). At the instability threshold, a dimensionless mea- 
sure of the ring self-gravity is GMd/if^n?') ~ 0.1, where 
f — 4 kpc and fl « 0.0458/to- The dependence of the 
growth rate is weU fitted by uj^to ~ 0.00842(Md- Md)"*^^, 
with the masses in units of 1O^"M0 and to = 10^ yr. For 

there is a "second instability" 



0.0196(Mrf - Mc2f-^^ 



Md > 3.12 X 10^° = Mc2, 

with growth rate Uito ~ O.vi-\)o{ivid — ivic2) 

The ratios of the complex perturbation amplitudes can 
readily be obtained from (62) - (65) once the 8, possibly 
complex frequencies are known. For the "first instability," 
for example, for Md — lO^^M© and the same conditions 
as for Figure 4, we find £1 k. 0.254 - 0.1 17i, Ai = 1 
(by choice), £2 « 0.0108- 0.365i, and A2 w 0.673+ 1.30i. 



These values correspond to </?i w 24.7°, tf2 ~ 91.7°, i/ii = 
(by choice), and V'2 ~ 297°. Thus the azimuthal density 
enhancement in the outer ring trails the density enhance- 
ment of the inner ring. On the other hand, the radial shift 
of the outer ring leads the shift of the inner ring. 

As a second example, for Md = 4 x 10"'^''Mq both the 
"first" and "second" instability occur. For the "first" in- 
stability, we again find that the azimuthal density enhance- 
ment of the outer ring trails that of the inner ring, whereas 
the radial shift of the outer ring leads that of the inner 
ring. For the "second" instability the situation is differ- 
ent in that both the azimuthal density enhancement and 
the radial shift of the outer ring trail those in the inner 
ring. Specifically, we find ipi « 178°, ip2 ~ 89.3°, V'l = 
(by choice), and ip2 ~ 283°. Note that for both rings, the 
angle of the radial shift is roughly 180° displaced from 
the azimuthal density enhancement. The displacement 
of the center of mass of the ring (equation 20) is domi- 
nated by the azimuthal displacement {\£i\ « 0.70|Ai| and 
\£2\ « O.67IA2I « 0.22|Ai|). Thus having ip^ and V'j about 
180° out of phase allows the center of mass of each ring to 
be closer to the origin, and this is a lower energy configu- 
ration. 

Consider now the angular momentum of the perturbed 
two ring system which from (58) is Vc/, — Vi^ + 1^24, = 
const. For the evaluation of P^, note that for each of the 
eight modes, Aj = —iuOaAj, and £j — —iuja£j, where 
q; = 1, .., 8 labels the mode. This implies that 



Wa = V'j + i{6]/5j) == Vj + i{<^]hj) 



(66) 



Thus, for the "first instability," for Md = W^Mq where 
ujto ~ 0.0500 + 0.00719i, we have V'j = <^j = 0.05/io, 
(which corresponds to forward precession), and 5j/5j = 
ij/ej ~ 0.00719/io- For an unstable mode, the coefficients 
of the six terms in Vip are all grow exponentially. The only 
possible way in which P^ =const can be maintained is to 
have V^ = = P^i + P^2- This relation provides a useful 
check on the correctness of the calculations. 

For the "first instability" {Md > Md), we find by eval- 
uating (58) that 7^01 — ~-V(f,2 > 0. This means that 
the angular momentum of the inner ring increases while 
that of the outer ring decreases. Thus the instability 
acts to transfer angular momentum inward. In contrast, 
for the "second instability" {Md > Mc2), we find that 
Vcf,! — —V,f,2 < 0. Thus, the "second instability" acts to 
transfer angular momentum outward. If P^i decreases and 
P^2 increases, then the average radius of ring 1 must de- 
crease and that of ring 2 must increase. Therefore, the 
"second instability" may be important for the accretion of 
matter in a gravitating disk. 

5. ECCENTRIC MOTION OF ONE RING AND "CENTER" 

Here, we consider the eccentric motion of a "disk" of 
one ring including the influence of the eccentric motion of 
the "center" Mq which is located at r = in equilibrium. 
The mass Mq includes the mass of a central black hole 
Mhfi if it is present. The ring perturbation is described by 
£{t) and A{t) as given by (52) and (53), while the "cen- 
ter" is described by £o{t) which is given by (49). With the 
perturbations ex exp(— iwt), we find 

[{lu -n)^ + E]£ + [2n{Lu -n) + d]A^ 2nl£o , 
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Fig. 4. — The figure shows the eccentric instability of a "disk" of mass Md consisting of two rings of equal mass Md/2 at 
equilibrium radii ri = 3 and r2 = 5 kpc. The top panel (a) shows the dependence of the real parts of the frequencies w^a 
(a = 1, .., 8) on Md- The labels ul, u2 indicate unstable branches where the frequency is complex; ul is referred to in the 
text as the "first instability" and u2 the "second instability." Here, to = 10® yr. The bulge has Mi, = 0.5 x IO^^Mq and 
rb = 1 kpc, and the halo Vh = 250 km/s and r^ — 5 kpc in the expressions given in §2.1. The bottom panel (b) shows 
the dependence of the growth rate uji on the ring mass. The onset of instability corresponds to the merging of the two 
real frequencies in panel (a). The tidal coefficients, {Cjk},, etc., are obtained using the equations of the Appendix and 
Avj = 2 kpc. 



[2n{uj -n) + d]£ + [{uj-nf + D]A 



^l£o 



(67) 



where fll = GMo/rf and f7g = GMi/rf, with Mi and ri 
the mass and radius of the ring. For a non-zero solution, 
the determinant of the 3x3 matrix multiplying (f , A,£o) 
must be zero. This leads to a sixth order polynomial in lo 
or w = tu — fl which can readily be solved (with Maple, 
R5) for the frequencies of the 6 modes oja (oi = 1..6). We 
obtain 



{w + nf -f7o^ 



{w^ 



E){r 



D) 



{2nw + dy 






5w^ 



E 
2D + - 



2Qw 







(68) 



where the strength of the interaction between the ring and 
the "center" is measured by 



ni 



flnflh 



Gy/MoMi 



(69) 



Here, D, E, d and il = ili are defined in (59), and ilg is 
defined in (51). 

Figure 5 shows the dependence of the growth rate uji = 
Tm{uja) on the mass of the "center" Mg with the mass of 
the ring held fixed. The associated real part of the fre- 
quency is positive. The onset of instability corresponds 



to the point where two of the six modes with frequencies 
given by (68) merge. The merging is again of positive and 
negative energy modes. The growth rate has to a good ap- 
proximation the dependence Wito ~ 0.00764(Mo--2.61)^/^, 
with Mo in units of IO^Mq and to = 10'^ yr. 

The ratios of the complex amplitudes follow from (67), 
and for the unstable mode for Mq = 2.75 x IO^Mq we 
find ujto « 0.155 + 0.00289i, nto w 0.131, ^oto « 0.142, 
nab ~ 0.0465, £ w -0.817 + 0.0333i, A == 1 (by choice), 
and £o « -3.38 + 0.916z which give \£\ w 0.817, \£i\ « 
3.51, ip w 182°, and (po « 195°, where to = 10*^ yr. Thus, 
the radial shift of the ring is roughly 180° away from the 
maximum of the density enhancement which has ip — 
(by choice). The shift of the center of mass of the ring 
(equation 20) is dominated by the azimuthal density en- 
hancement. Thus the radial shift and azimuthal displace- 
ments are such that the center of mass of the ring moves 
closer to the origin which is a lower energy configuration. 
Note that the radial shift of the "center" trails the ring 
center of mass by an angle 360° — ipo ^ 165°. Hence, the 
torque of the "center" on the ring acts to reduce the angu- 
lar momentum of the ring as verified below. At the same 
time the ring acts to increase the angular momentum of 
the "center." 

Consider now the angular momentum of the perturbed 
ring plus "center" system which is given by (58). Fol- 
lowing the arguments of the prior section, the sum of the 
angular momentum of the "center" and that of the ring 
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Fig. 5. — The figure shows the eccentric instability of a ring of mass Mi interacting gravitationally with a displaced 
point mass Mg shifted from its equilibrium position r = 0. The halo and bulge potentials are the same as in Figure 1. 
The mass Mq includes the mass of a central black hole Mbh if present. The mass of the ring is Mi = 6 x 10^'^Mq and its 



radius is ri 
Ari = 2 kpc. 



3 kpc. Here, to = 10 yr. The tidal coefficients are evaluated using the expressions of the Appendix with 



must be zero for a growing mode. The angular momen- 
tum of the "center" is simply MgeQipo- For a pure mode, 
we have TZe{uja) ~ ifo- As mentioned, the real part of the 
frequency is positive for the unstable mode and therefore 
the angular momentum of the "center" increases while the 
angular momentum of the ring decreases. Thus, there is a 
transfer of angular momentum from the ring to the "cen- 
ter." Due to the loss of angular momentum the average 
radius of the ring will decrease. Thus the instability may 
be important for accretion of matter to the galaxy center. 



to give ipj (t) , which is the angle to the maximum of the 
radial shift, and tpjit): which is the angle to the maximum 
of the azimuthal density enhancement. These angles are 
analogous to the line-of-nodes angles for the tilting of the 
rings of a disk galaxy (Lovelace 1998). Thus, we solve 
lOiV first order equations. In all cases, the total energy 
(54) and total canonical angular momentum (58) are ac- 
curately conserved. The different frequencies Qj, ilgj, and 
flAj, and the tidal coefficients {Cjk}, etc., are evaluated 
using the equations of the Appendix. 



6. ECCENTRIC MOTION OF N RINGS 

For the results presented here, the rings are taken to be 
uniformly spaced in r with radii rj = 1 + {j — l)Sr kpc 
with 6r = 0.5 kpc and with j = 1, .., A^ = 31. The value 
A" = 31 gives good spatial resolution over all but the inner 
part of the disk. The outer radius rjv (in the range say 
10 — 20 kpc) has little influence on the eccentric motion de- 
scribed here, as verified by comparing results with r^ = 16 
kpc those obtained with significantly larger r^. Also, the 
eccentric motion of the outer disk, say, r ^ 3 kpc, is es- 
sentially independent of dr. We first consider in §6.1 the 
case where the ring masses correspond to the exponential 
distribution discussed in §2.1. The inner part of the disk 
is found to be strongly unstable to eccentric motions and 
therefore in §6.2 we consider a disk with the mass of the 
innermost three rings reduced. In §6.3 we consider disks 
with a smooth reduction in T,d{r) in the inner part of the 
disk, r <^rd. 

We solve (52) and (53) numerically as eight first order 
equations for e^j, ixj, £yj, and iyj, and for Sxj, S^j, Syj, 
and 6yj, j = 1, .., N . At the same time, we solve the two 
additional equations. 



(70) 



6.1. Exponential Disk 
Here, we consider the eccentric motion of the rings 



dipj 


^xj ^yj ^yj ^xj 


dt 


e2 . + e2 . 


dxpj 


5xj5yj - SyjSxj 


dt 


5% + ^% 



for the case where Mo 



2'Krj5rY,d{rj) with Sd(r) 



(71) 



I](ioexp(— r/rrf). The mass of the center is assumed given 
by (47), which gives Mo « 1.06 x 10^ for the parameters 
of Figure 1. Alternatively, this value of Mq could be due 
in part to a central black hole. 

Figure 6 shows the dependences of the radial shifts 
Sjifj) and azimuthal displacements Sji^ipj) of the rings 
( j = 1 — 31) and the radial shift of the "center" eoi'Po) at 
a short time, 100 Myr after an initial perturbation. This 
type of plot is related to the plots emphasized by Briggs 
(1990) for characterizing the warps of galactic disks (see 
also Lovelace 1998). The angles (p and V" are analogs of 
the line-of-nodes angle the warp. 

From Figure 6 note that the azimuthal displacements Sj 
are larger than the radial displacements ej so that the dis- 
placement of the center of mass of a ring is dominated by 
Sj. The eccentric motion of the inner rings, say, j = 1 — 4 
or ri = 1 to r4 = 2.5 kpc, show the most rapid, expo- 
nential growth. For these rings the angles ipj and tpj are 
approximately 180° degrees out of phase, and this agrees 
with the behavior found for the "second" instability of two 
rings discussed in §4. As mentioned, this allows the center 
of mass of each ring to move closer to the origin, which 
is a lower energy configuration. Note that with increasing 
j, both ipj and ipj decrease for the inner rings which cor- 
responds to a trailing pattern the same as found for the 
second instability of two rings. Note also that the shift 
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Fig. 6. — Polar plot of the radial shift ej and azimuthal displacement Sj of ring matter as a function of the angles (pj 
and V'j (j = 1 — 31) at time t = 100 Myr. The radial shift of the "center" eq is indicated by the solid dot. The conditions 
correspond to the galaxy parameters of Figure 1. The rings have radii Vj = I + Sr{j — 1) kpc and 5r = 0.5 kpc for 
j = 1,..,31, masses Mj = 2'Krj5rYjd{rj) with E^ given in §2.1, and Mq given by (47), which gives Mq « 1.06 x IO^Mq. 
The initial values of the shifts and displacements are e^ = 0.1{rj/rmax) = ^j, fj=^ = ''Pjj ^^'^ ^Ox = 10^^ and egy = 0. 
The units of ej and Sj are arbitrary in that the equations are linear. 



of the "center" eo trails the shift of the center of mass of 
the j = 1 ring in agreement with §5. Thus the torque of 
the first ring on the center acts to increase the angular 
momentum of the center while the torque of the center on 
the first ring decreases the rings angular momentum. 

Figure 7 shows the exponential growth of the azimuthal 
displacement Si and radial shift ei of the first ring and 
the simultaneous growth of the radial shift of the center 
eo- The e— folding time is about 29 Myr. For comparison, 
the period of oscillation of the center is Tq = 27r/J7o ~ 46 
Myr for the conditions shown, where f2o is given by (51). 
The growth of the eccentric motion of the inner rings is 
reduced somewhat if the mass of the center is reduced to 
Mq = 10^ Mq; the e— folding time for in this case is about 
38 Myr for ring 1. Figure 8 shows the perturbations of 
the angular momentum of the center Pq and the rings Pj 
at i = 100 Myr obtained from (58). In agreement with §5 



and the abovementioned direction of the torque, the an- 
gular momentum of the center increases while that of the 
first and second ring decrease. The decrease in angular 
momenta of these rings will result in their radii shrinking. 
Note that the center rotates in the same direction as the 
disk matter. 

The present linear theory does not address the issue of 
saturation of growth of the eccentric motion. One possi- 
bility is that the strong instability of the inner rings of the 
disk leads to the destruction of this part of the disk. 

6.2. Exponential Disk with Rings 1-3 Reduced 

Here, we consider the eccentric motion of the rings for 
the case where the ring masses are the same as in §6.1 ex- 
cept that Ml -^ lO-^Mi, M2 -^ O.IM2, and M3 -^ O.3M3, 
which are the rings with radii ri = 1, r2 = 1.5, and 
rs = 2 kpc. The mass of the center is the same as in 



Fig. 7. — The plot shows the exponential growth of center shift eq a-nd radial shift ei and azimuthal displacement Si of 
the first ring at ri = 1 kpc for the same conditions as for Figure 6. 
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Fig. 8. — Plot of the perturbations of the angular momentum of the center Pg ^nd the rings Pj at t = 100 Myr for the 
same conditions as Figure 6. 



§6.1, Mo ~ 1.06 X IG^M©. The disk mass is reduced by 
a factor « 0.84 compared with an exponential disk. The 
aim of reducing Mi — M3 is to reduce the growth rate of 
the eccentric motion of this part of the disk. 

Figure 9 shows the essential behavior in a polar plot of 
the displacements and shifts at two times. The radial shift 
of the center is negligible on the scale of this figure. The 
e— folding time for ring 3 is about 49 Myr. Notice that 
the curves traced out by Sj and by ej are approximated 
straight lines from the origin which rotate rigidly in the di- 
rection of motion of the matter for j = 4 to about j — 11, 
which corresponds to r^, = 2.5 to rn = 6 kpc. The instan- 
taneous period of rotation of this pattern is « 60 Myr, 
which is longer that the oscillation period at the center, 
Tq w 46 Myr. This case is an example of the phase-locking 
of the eccentric motion of these rings due to the self-gravity 
between the rings. This phase-locking is analogous to that 
which occurs in the tilting motion of the rings represent- 
ing a disk galaxy due to self-gravity (Lovelace 1998). In 
the case of tilting of rings the phase-locking results in the 
line-of-node angles of the rings in the inner part of the disk 
becoming the same. 

Also in this case the growth of the eccentric instability of 
the inner rings is sufficiently fast that it probably further 
disrupts the inner part of the disk. 

6.3. Reduced Inner Disk 

Here, we study the eccentric motion of a disk the in- 
ner part of which is attenuated relative to an exponen- 
tial disk. Specifically, the rings masses are Mj — 2TTrj6r 
Sd(r,), where 



T.d{r) = Sdocxp (-13—) exp 

r / \ rd 



(72) 



with /3 =const. We consider j3 = 1. The mass of the center 
is the same as in §6.1, Mq « 1.06 x IO'^Mq. However, the 
motion of the center is negligible and value Mq has little 
influence on the eccentric motion of the disk described be- 
low. The disk mass is smaller than for an exponential disk 
by a factor w 0.47. 

Figure 10 shows the essential behavior in a polar plot 
of the shifts and displacements at two times. The shift 
of the center is negligible on the scale of the plot. The 



magnitude of azimuthal displacement of the outer ring 
5m exhibits an approximately linear growth with time, 
5m ~ const -I- i/660Myr, for the considered initial con- 
ditions and t ^ 1 Gyr. In contrast, the magnitude of the 
radial shifts ej remains bounded by its initial maximum 
value eAr(t = 0). 

The patterns formed by both 5j and ej in Figure 10 
are trailing spirals. This is different from the proposal of 
Baldwin et al. (1980) that leading spirals should form. 
The rotation of the outer point on the spiral Ajv(rAr = 16 
kpc) in Figure 10 is in the direction of rotation of the 
disk matter. Its pattern speed Jlp corresponds to a period 
2tt /Vtp K, 460 Myr, which is a factor k, 1.24 longer than 
the rotation period of matter at this radius (w 370 Myr). 
The pattern period of say Ai5(r = 8kpc) is ss 188 Myr, 
which is a less than the rotation period of the matter at 
this radius (« 206 Myr). Thus, it is evident that the spiral 
is "wrapping up" as time increases. At the same time, the 
spiral pattern propagates radially outward. The outward 
speed is about 10 km/s at r ^ 8 kpc for t ^ 300 Myr. 

Figure 1 1 shows a polar plot of the radius to the max- 
imum of the azimuthal displacement rj{tpj) at two times. 
The curve is a trailing spiral with an approximate fit given 

by 



if) = A{t) exp 



a{t) 



(73) 



where A w 0.065 (i/kpc) rad, and a « 7.0+0.0044 (i/Myr) 
kpc for i <^ 1 Gyr. Thus the radial spacing between spiral 
arms is Ar ~ {2'iTa/A) exp(r/a) for A^ ^ a. For validity of 
the ring representation we must have Ar > 25r{= 1 kpc), 
where 5r the separation between rings. 

Figure 12 shows the radial variation of the perturbations 
in ring angular momenta. The perturbation of the angular 
momentum of the central mass is negligible. This figure 
should be compared with Figure 8 which gives the same 
plot for an exponential disk. 

Some simplification of equations is possible in the 
present case at long times due to the fact that 5j 3> ej. 
Firstly, we have (5v k, 5va,4> with 



5v^{r, 4>) w —{5x + fl5y) sin0 + {dy — il5x) coscp . (74) 
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Secondly, 



E(r) 



- (Sx cos (/) 



Sy sin ( 



(75) 



for 6j » rcj/rd- 

Figure 13 shows the profiles along the x— axis through 
the galaxy center of the fractional change in the surface 
density S'E/'E and the change in the aziniuthal velocity 
Sv^ obtained from (74) and (75). The opposite signs of 
6v^ on the two sides of the galaxy would of course make 
the rotation curves on the two sides different as observed 
in some cases (Swaters et al. 1998). Note that in some re- 
gions the changes (JS and 5v^ are correlated and in other 
regions they are anticorrelated. Figure 14 shows two- 
dimensional appearance of the fractional surface density 
variations from (75). 

For long times t > \ Gyr, the the azimuthal displace- 
ments and shifts of the inner rings (2 and 3) start to be- 
come large compared with the values in the outer disk 
(r > 4 kpc) even though these rings have very small 
masses. At the same time, the displacement of the cen- 
ter, which has mass Mq = 1.06 x lO^M©, grows, and at 
t = 1 Gyr it is eo ~ 0.057 for the conditions of Figures 
9-13. If the mass of the center is Mq = lO^M©, then the 
displacements and shifts of rings 2 and 3 at t = 1 Gyr are 
significantly reduced as is the shift of the center which is 
eo « 0.014. 

7. CONCLUSIONS 

The paper develops a theory of eccentric {ma = ±1) lin- 
earized perturbations of an axisymmetric disk galaxy re- 
siding in a spherical dark matter halo and with a spherical 
bulge component. The disk is represented by a large but 



finite number N of rings with shifted centers and with per- 
turbed azimuthal matter distributions. This description is 
appropriate for a disk with small 'thermal' velocity spread 
vth where the matter is in approximately laminar circular 
motion. The spread for a thin disk has (wt/i/w^)^ <C 1, but 
it is sufRcent to give a Toomre Q{r) ^ 1. Earlier, Bald- 
win et al. (1980) discussed asymmetries in disk galaxies 
in terms of shifted rings but without interactions between 
the rings and without azimuthal displacements of the ring 
matter. Account is taken of the shift of the matter at the 
galaxy's center, which may include a massive black hole. 
The gravitational interactions between the rings and be- 
tween the rings and the center is fully accounted for, but 
the halo and bulge components are treated as passive grav- 
itational field sources. Equations of motion are derived 
for the ring and the center, and from these we obtain the 
Lagrangian for the rings+ccnter system. For this system 
we derive an energy constant of the motion, and a total 
canonical angular momentum constant of the motion. 

We first discuss the nature of the precession of a single 
ring with the other rings fixed; this case although not self- 
consistent is informative. There are four modes, analogs 
to the normal modes of a non-rotating system, and two 
have negative energy and two positive energy. Negative 
energy modes are unstable in the presence of dissipation 
such as that due to dynamical friction. We go on to study 
the eccentric motion of a disk consisting of two rings of 
different radii but equal mass Md/2. Above a threshold 
value of Md the two rings are unstable with instability 
due merging of positive and negative energy modes. This 
result is obtained by solving the eighth order polynomial 
for the frequencies of the eight modes. Above a second, 
somewhat larger threshold value of M^, a second instabil- 
ity appears, and in this case the ring motion is such that 
the angular momentum of the inner ring decreases while 



Fig. 9. — Polar plot of the radial shift ej and azimuthal displacement 5j of ring matter as a function of the angles ipj 
and V'i of the maxima of the shift and displacement at times t = 300 and 305 Myr. The solid dot labeled by the arrow 
indicates the shift of the center at t = 300 Myr; the other dot is the shift at 305 Myr. The rings, the center, and the 
initial values of the shifts and displacements are the same as in Figure 6 except that Mi -^ 10~^Mi, M2 -^ O.O3M2, and 
M3 -^ O.IM3. Thus the conditions correspond to the galaxy parameters of Figure 1 except that the disk mass is reduced 
to « 5.06 X 10^'^Mq. The shifts and displacements of the rings j = 1, 2 are dynamically unimportant and are not shown. 
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Fig. 10. — Polar plot of the radial shift Cj and azimuthal displacement 5j of ring matter as a function of the angles 
tpj and V'j of the maxima of the shift and displacement at times t = 300 and 400 Myr. The radial shift of the center is 
negligible on the scale of this plot. The rings, the center, and the initial values of the shifts and displacements are the 
same as in Figure 6 except that the ring masses are obtained from (72). The mass of the center is the same as in Figure 
6, Mg ss 1.06 X IO^Mq. Thus the conditions correspond to the galaxy parameters of Figure 1 except that the disk mass 
is reduced to w 2.82 x IQ^^Mq. The shifts and displacements of the rings j = 1, 2 are dynamically unimportant and are 
not shown. 



that of the outer ring increases. For the unstable motion, 
the maximum of the azimuthal density enhancement of a 
ring occurs at an angle about 180° from the direction of 
the radial shift. This allows the center of mass of the ring 
to move closer to the center of mass of the other ring and 
to the origin. 

We also analyze the eccentric motion of a disk of one 
ring interacting with a radially shifted central mass. This 
system has six modes, the frequencies of which are ob- 
tained by solving a sixth order polynomial. In this case, 



instability sets in above a threshold value of the central 
mass (for a fixed ring mass), and it acts to increase the 
angular momentum of the central mass (which therefore 
rotates in the direction of the disk matter), while decreas- 
ing the angular momentum of the ring. The instability is 
again due to the merging of positive and negative energy 
modes. 

We study the eccentric dynamics of a disk with an expo- 
nential surface density distribution represented by a large 
number TV = 31 of rings and a central mass Mq ~ IQ'^Mq 



Fig. 11. — Polar plot of the radius to the maximum of the azimuthal displacement rj{'ipj) at two times for the same case 
as Figure 10. 
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Fig. 12. — Radial dependence of the perturbations in the ring angular momenta Pj at t = 300 Myr from equation (58). 
The angular momentum of the center Pq is negligible. 



which may include the mass of a black hole. The outer ra- 
dius of the disk is r^ — 16 kpc; we have checked that this 
value has negligible affect on the reported results. In this 
case, we numerically integrate the equations of motion. A 
check on the validity of the integrations is provided by 
monitoring the mentioned total energy and total canoni- 
cal angular momentum, which are found to be accurately 
constant in all presented results. The inner part of the 
disk r ^ 2.5 kpc is found to be strongly unstable with 
e— folding time ^ 30 Myr for the conditions considered. 
The e— folding time is somewhat longer if Mq = 0. An- 
gular momentum of the rings is transferred outward, and 
to the central mass if it is present. A trailing one-armed 
spiral wave is formed in the disk. This differs from the 
prediction of Baldwin et al. (1980) of a leading one-armed 
spiral. The outer part of the disk r ^ r^ is stable and in 
this region the angular momentum is transported by the 
wave. Thus our results appear compatible with the theo- 
rem of Goldreich and Nicholson (1989) regarding angular 
momentum in stable rotating fluids. The instability found 
here appears qualitatively similar to that found by Taga 
and lye (1998b) for a fluid Kuzmin disk with surface den- 
sity S ex 1/(1 + r^)'^/^ with a point mass at the center 



where unstable trailing one-armed spiral waves are found. 

The present linear theory does not address the issue of 
saturation of growth of the eccentric motion. One possi- 
bility is that the strong instability of the inner rings of the 
disk leads to the destruction of this part of the disk. For 
this reason we have studied a disk with a modified expo- 
nential density distribution where the surface density of 
the inner part of the disk is reduced. However, the mass 
of the center of the galaxy was kept the same as in the 
case of an exponential disk, Mq ~ lO^M©. In this case 
we find much slower, linear - as opposed to exponential - 
growth of the eccentric motion of the disk for times t <^ 1 
Gyr. A trailing one-armed spiral wave forms in the disk 
and becomes more tightly wrapped as time increases. An- 
gular momentum is transferred outward. The motion of 
the central mass if present is small compared with that of 
the disk for t ^ 1 Gyr. 

For long times t > 1 Gyr, the the azimuthal displace- 
ments and shifts of the inner rings start to become large 
compared with the values in the outer disk. At the same 
time, the radial shift of the center grows. This shift is 
significantly reduced if the mass of the center is changed 
from - 10^' to IQ^Mq. 



Fig. 13. — Profiles of fractional density variations 10^(5S(a;, 0)/S(a;, 0) (in percent) and variation of azimuthal velocity 
lQ^6vy{x,())/v^{N) along the x— axis through the middle of the galaxy at t = 400 Myr for the same case as Figure 10. 
Here, v^{N) « 265 km/s is the disk rotation velocity at r = 16 kpc. The vertical scale is arbitrary in that the equations 
solved are linear, but the ratio S'E/Svcf, is fixed. 
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Fig. 14. — Two-dimensional appearance of the fractional density variations \Q'^5Yi{x^y)/Yi{x,y) at i = 400 Myr for the 
same case as Figure 13. 
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APPENDIX 
TIDAL COEFFICIENTS 



For \rj — rk\ ^ -v/Arj Ar^, the ring profiles can be treated as delta functions, S{r\rj) -^ 5{r — rj)/r, and consequently 
the 'tidal coefficients' of equations (29) - (32) can be simplified to give 



Cjk ^ 


a GMjMk 


D,, . 


a GMjMk 


D',, ^ 


a GMjMk 


Ejk ? 


a GM.Mk 



drj dvk 
dlC(rj,rk) _ 
Tk drj 

dlC{rj,rk) 

rj dvk 
fC{rj,rk) 



2GM,-Mfe 



GMjMk 



E{kjk) , 



7rr2(r2-r2) 



= Dkj , 

GMjMk 



{tj + rk)E{kjk) + (rj - rk)K{kjk) 



{r]+rl)K{k,k)-{rj+rkfE{k,k) 



Tj rk '^r'^rKrj + Vk) 

(Mathematica V. 3) where kjk ^ 2^/rjrk/{rj + r^), where JC is defined by equation (31), and where 

r.7r/2 / I i.tt/2 



(Al) 
(A2) 
(A3) 
(A4) 



K{k) 



b/^Jl-k^sm^<j> , and E{k) = f d<P\J I - k^ sin^ , 



are complete elliptic integrals of the first and second kinds respectively. Note that for r^ ^ r^, Cjk ~ —GMjMk/r^, 
Djk « GM,Mk/(2rl), D'^^ « -GM,Mk/rl and E,k « GMjMk/(2rl). 

In the opposite limit where \rj — r^j <C ^tj rk , equations (27)-(30) can be evaluated approximately as 

GM.Mk r ydy exp(-yV4) 



a 



jk 



87r-/7r(Ar)r2 7_^ rk — r-j + yAr 



2f — (j'A; — ''j + J/Ar) In 



kfe - rj + yAr 



GMjMk 
27r(Ar)2f 



1 - ^uexp(-uV4)crfi(u/2) + ©(Ar/f) 



(A5) 



where r = (r^ + rfe)/2, Ar = W (Ar| + Ar^)/2, u = (r^ — rj)/Ar, and and erfi(x) = crf(«a;)/« = (2/^/71) J^ dyexp{y'^). 

The integral in (A5) is a principal value integral of the form occurring in the plasma dispersion function W (Ichimaru 
1973). Also, 

dy cxp(— y^/4) 



D 



jk 



GMjMk 
GMjMk 



rk — rj + yAr 

dy exp(— 1/^/4) In 



2r — {rk — rj + yAr) In 
8f/Ar 



8(f/Ar) 



and 



E. 



jk 



GMjMk 
2'K\rKr^ 



dy exp{—y /4) In 



u + y\ 
- 27r(f/Ar) exp(-MV4)erfi(u/2) 



1.0827(f/Ar) 



\u + y\ 



From these expressions we obtain 



G 



GMf 



n 



2Trrj{Arj 



D 



n 



GM] 
27rr3 



In 



10.68rj 
Ar,- 



E. 



. GM] 



■In 



1.445rj 

Ar, 



(A6) 
(A7) 

(A8) 



Equations (Al) - (A8) are valuable for numerical evaluation of the tidal coefficients. 

In the following we derive some useful relations involving the tidal coefficients. From equation (12) we have 



^\{r) = i^ = 2^G|"r'dr'Erf(r')^ 
The ^ integral in this case is 



d* 



(1 - r'cos*/r) 



27r fT-2 



j,2 ^ {^^^Y — 2rr' cos ^ 



3/2 



Comparison with (A2) shows that shows that 



— Ur + r')E+{r-r')K] 



^% 



rdrS{r\rj)nl{r) 



J fc=i 



(A9) 
(AlO) 

(All) 
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This expression does not include the disk mass within the inner ring. As discussed in §2.6, this part of the disk is treated 
as a point mass Mq with unperturbed position r = 0. If there is a central black hole M^^, its mass is include in Mq. To 
account for the influence of Mq, we simply add the term GMo/r^ to the right hand side of (All). The resulting expression 
for Udj is useful for the numerical calculations. For TV ^ 30 and rj = 1 to 10 — 20 kpc, we find that equation (All) gives 
accurate agreement with the analytic expression (3). 

An alternative expression for fl^ can be obtained by integration by parts, 



^lir) = - 



2ttG 



/ r'dr'^ IC{r,r') 
Jo or' 



Thus 



or 



M,. I M 

or 



rdrS{r\rj)r 



dnjir) 



GMj ^Mk I I rdr r'dr'S{r\rj) 



dS(r'\rk)d[]C{ry)/r] 



dr 
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^ ) = GMj ^Mk j j rdr r'dr'S{r\rj) 



1 d[r'S{r'\rk)] S{r'\rk) 



J k 

In view of equations (27) - (30), this equation can be written as 



dr' 



IC{r,r') dIC{r,r') 



dr 



M, I rfl 

or 



= Y,{C,k + D,k) - Y^^D'^k + E,k) 



This relation is useful in §2.5. 

We also evaluate the disk gravitational potential in the ring representation, 
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-G / d\' Ml) 



r — r 
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where 



$ 



dj 



7r(r + r') 

o 

rdr S{r\rj)^d{i') 



wX""'' 



For \rj - rk\ > ^/Ar~Ar^, 
whereas for \rj — rk\ ^ y/i~ji^k 



Ajk = -GM.Mk / / rdr r'dr' 5(r|r,)5(r'|rfc) 
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(A19) 



where u= {r^ — rj)/Ar as above. Note that Aj^ = A^j and that Ajj w —GM"^ {1 / nrj) ln(lQ.68rj / Arj) . 
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